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Abstract 

We investigate the dependence of s-wave meson-baryon scattering amplitudes on different regular- 
izations within the framework of the chiral unitary model. We employ two different regularization 
schemes, i.e. dimensional and form-factor regulariaztions to tame the divergences in the model. 
We also study the analytic structures of T-matrices, using those regularization schemes. We find 
that while the form-factor regularizaion produces almost the same results as the dimensional reg- 
ularization did, the on-shell approximation is to some extent limited in the case of the form-factor 
regularization. Having chosen parameters properly, we show that the regularization dependences 
can be minimized. 
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I. INTRODUCTION 



Understanding meson-baryon scattering has been a very important issue for several 
decades, since it gives information not only on the strong interaction between hadrons but 
also on the origin of baryonic resonances. Recently, chiral perturbation theory has been of 
great success in explaining low-energy meson-baryon scattering, in particular, strangeness 
S = 0, — 1 channels. While 5- wave vrA^ and K^N scattering can be well described by the 
Lagrangian of leading order in the chiral expansion, the KN system requires multi-coupled 
channels in order to generate resonances such as A(1405) P,!^!- Though chiral perturbation 
theory (ChPT) provides us with a theoretical and systematic framework to study meson- 
baryon systems, it is restricted to lower energy regimes. In order to explain the KN system 
quantitatively, one has to go beyond ChPT. 

Recently, there have been noticeable works on the KN systems, based on the effective 
chiral Lagrangian: Kaiser et al. P| examined the S'-wave KN system. They utilized three 
different types of the pseudo-potentials, keeping the low energy constants from the effective 
chiral Lagrangian and solved the Lippmann-Schwinger (LS) equation in the coupled- channel 
formalism. Krippa and Londergen jj] investigated also KN scattering, using the on-shell 
approximation. Though Ref. seems to describe the experimental data 0, |J well, thresh- 
old behavior can not be reproduced well in the S^vr^, S+tt^, and Avr^ channels due to 
the openings of new channels. The Valencia group also studied extensively the s-wave 
system Q. Ref. 3 started from the effective chiral Lagrangian to construct the pseudo- 
potentials and used the on-shell factorization to solve the coupled-channel LS equation, 

n 

employing a regularization with the three-dimensional cut-off. Later, [8] introduced the di- 
mensional regularization instead of the three dimensional cut-off. They also made use of 
the N/D method to solve the scattering equation. Lutz and Kolomeitsev Q investigated 
meson-baryon scattering in the context of the effective chiral Lagrangian with the large N^ 
counting. They introduced a minimal chiral subtraction scheme within the dimensional reg- 
ularization in order to keep the covariant chiral counting rules. Though there exist technical 
differences between these works, almost all these works showed a remarkable agreement with 
experimental data. Thus, it is of great significance to understand some essence, if any, in 
dynamics of different technical methods, in particular, different regularization schemes and 
the role of corresponding parameters. 
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In the present work, we turn our attention to the dependence of the s-wave KN and 
vrA^ systems on different regularization schemes. While the dimensional regularization is 
preferably selected in ChPT, the conventional meson-exchange model jl^ |^ gets used to 
introduce form factors, which describe the extended hadron structure, as a regularization. 
Though the dimensional regularization gains an advantage over the form factors in many 
aspects, in particular, in the context of renormalization and gauge invariance, it is less 
convenient to associate with a LS-type scattering equation. Hence, in this work, we will 
show that the use of form factors presents almost the same results as that of the dimensional 
regularization. The present investigation will shed light on the meaning of the parameters 
involved in the regularization and will pave the way for solving the LS scattering equation 
without on-shell approximation in the future work by employing the form factors in place 
of the dimensional regularization. 

The paper is organized as follows: In section II, we shall explain the formalism of the 
chiral unitary model to investigate the meson-baryon systems. The analytic structure of the 
scattering amplitudes will be presented. In section III, we shall give the numerical results 
in the S = channel with two different regularizations used. In section IV, we shall present 
those in the S = — 1 channel. In section V, we shall summarize the present work and draw 
conclusion with outlook. 



II. FORMALISM 



We start with the effective chiral Lagrangian to the lowest order C\[q [12 



where 



C\jB = {B(zf - M)B) + -D{BYl,{u„B}) 
+ ^F{BY75K,B]), 



(1) 
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B and $ represent baryon and pseudoscalar fields, respectively: 
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The coefficients F and D denote the reduced matrix elements for semiletonic decays of octet 
baryons in SU(3) flavor symmetry. / in the U field is the meson decay constant and 
(■ ■ ■ ) denotes the trace over SU (3) flavor space. Since we are interested in s-wave scattering 
at low energies, we need only the terms of order 0{p) in the expansion of the U field: 

The Lagrangian in Eq. (0} is also known as the venerable Weinberg- Tomozawa term. The 
corresponding pseudo-potentials can be easily obtained as: 



4/^ 



'3 '■> 



(5) 



where pi and pj {ki and kj) are initial and final baryon (meson) momenta, respectively. 
The subscripts i and j denote the indices representing the coupled-channel states. The 
coefficients Cj,- are derived from Eq. (HI). Explicit forms of the Ca for S = — 1 and S = are 
given in Refs. [7|, HJ] . 

The next step is to solve the Bethe-Salpeter(BS) equation. The general T matrix in the 
coupled-channel formalism is given by 

A Vu{A + Mi)Ti, 



(6) 



, . {2nr {{P - qY - Mniq' - mt J 
where M/ and mi represent the octet baryon and meson masses in the intermediate state, 
respectively. P designates the total momentum given hj p + q = {^/s, 0, 0, 0) in the center 
of mass system, where p and q denote the intermediate baryon and meson momenta, respec- 
tively. Having utilized the on-mass-shell factorization or N/D method, Eq. (jHl) becomes an 
analytically solvable and pure algebraic equation: 

d^q 2Mi 



T- 




{27rf {{P - qf - Mf}{q^ 




(7) 



The pseudo-potential in Eq. ((Tj) is then written by 



V(^), = - M. - A« ^^^^ (8) 

where Ei stands for the zth baryon energy. 

In order to solve the BS scattering equation, we have to introduce the regularization to 
remove the divergence. We first use the dimensional regularization (DIM) as in Refs. M, 

nn 

convenience without loss of generality, the indices representing the coupled 
channels will be omitted from now on. Employing the dimensional regularization, we obtain 
the familiar result for the two-body propagator q]: 

2M fm^-M^ + s, m2 ^ + - s - n 2M , 



where 



^ = y^(M2^^^^j2T^^2^r^^ ^ - (M - m)2) (s - (M + m] 



2) 



related to on-shell center of mass momentum of meson-baryon system, ^/(2a/s). Also it is 



related to the phase space (p = ^) or the flux factor which governs the imaginary part 
of the T-matrix. /i stands for the renormalization scale, which contains information on the 
divergence of G{^/s) with the minimal subtraction. 

Now, we introduce form factors instead of the dimensional regularization. The form 
factors encode the complicated extended structure of hadrons and are usually parameterized 
in the form of: 

= (10) 

where A denotes a four dimensional cut-off parameter. In the case of n = 1, i.e. monopole 
type of the form factor (MF), it is equivalent to the well-known Pauli-Villars regularization. 
The loop integral with the monopole type of the form factor can be performed analytically: 

^, - 111 '^J± (A^ - 

C.(V^)MF j _ _ ^2}(^2 _ ^2) (A2 - q2) 

2M f - + s , f + - s - ^ 

In — - H In 
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}, (11) 



5 



where 

\2 



r/ = y(M2-m2-s) -4sA2. (12) 

Though the G{y/s)MF shows a similar structure to G(v^)dim in Eq. 0, their analytic 
behaviors are rather different. It is straightforward to calculate all the cases of n > 1 by 
using the following recursion formula: 

G„+i(v/^) = G„(v^) - ^-^_G„(v/^). (13) 

For example, the analytic form of G'(v^)df for the dipole-type form factor (n = 2) becomes 
^. ^. ^, ^, Af(A2-m2) ri M2-m2-s + 2A2-r/\ 

G DF = G MF \-^''Tr2 2 ^oa2^ f" 

The divergence can be canceled again by introducing the counter terms with a subtraction 
parameter a: 



G^G+j^a, (15) 



2M 
16^ 

In coupled channel calculations, G and a are diagonal elements of a matrix in the meson- 
baryon channels. The subtraction parameter a in Eq. ()15|1 change only the real part of 
G{\^), while the imaginary part should be independent of regularization schemes, since it 
is finite as is constrained by unitarity. 

We plot the real part of G{^/s) with the DIM, DF, and MF for four different channels, i.e. 
irN and vrS without the subtraction parameters in Fig. ^ Here, we use /i = 1200 MeV for 
the DIM, whereas A = 1000 MeV for the DF and MF. Thus, we can minimize the difference 
the between regularization schemes by fitting the subtraction parameters. However, the 
results turn out to be different above the thresholds. 

In order to find resonances in each channel, we have to look for the corresponding poles 
in the second Riemann sheet. Since the pseudo-potential is real as given in Eq. (jSJ, the 
poles in the second Riemann sheet are effected only by G{^/s). The analytic structure of 
the two-body propagator is as follows: 

G{z)2nd = G{z)ist, for < M + m, 

G(2)2„d = G{z),,, - 2ilmG{z)i,t, for > M + m, (16) 

where z stands for the total energy on the complex plane. Since the cut starts from the 
threshold points, the imaginary part of G{z) can be obtained by the discontinuity along the 



cut y/s > M + m: 



ImG'(v^) 



P Mi 



(17) 



2 

H' he s 
. However, we have two complex variables, i.e. 



In the case of the DIM, we find exactl y t he same expression for the ImG'(y^) by taking 



^ and 7] give two different branch cuts which start from the points of ^/s = M + m and 
where Eq. p2j) equals zero in the region of 1 GeV ~ 2 GeV in the case of the form-factor 
regularizations. Hence, we have to take into account these two variables to find the poles 
in the second Riemann sheet. The variable 1] must be pure imaginary so that Eq. p7|) may 
be satisfied. While the on-mass-shell approximation given in Eq. ((7j) works perfectly well 
for the DIM, it is to some extent limited for the form-factor regularization. However, since 
we are interested in resonances in the region of 1 GeV ~ 2 GeV, in which t] becomes pure 
imaginary as shown in Fig. |21 the unitarity condition is also well satisfied in the present 
case. Therefore, the term with rj [t] term) in Eq. pi|) becomes pure real in the regiqri of 
1 GeV ~ 2 GeV. For experimental data of S* = and S = —1 we use Refs. jl7 
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III. NUMERICAL RESULTS OF S = MESON-BARYON SECTOR 

We start with the S = meson-baryon sector from which A^*(1535) (/ = 1/2) [2^ and 
A(1620) (/ = 3/2) resonances arise. In Table|l]we list the possible coupled channels. We fix 
the parameters to describe the S = sector as follows: We choose the meson decay constants 
according to the empirical data = 93.0 MeV, fr, = 120.9 MeV, and fx = 113.46 MeV. 
The renormalization point for the DIM is set tohe fi = 1200 MeV and the cut-off parameter 
is fixed to be A = 1000 MeV for the DF and MF. We use the subtraction parameters for 
the DIM, DF and MF as given in Table ITTl While we take the subtraction parameters for 
the DIM from Ref. |l^, we determine them for the MF and DF in such a way that the 
threshold-point values of G{^/s) are equal to those of the DIM in order to keep consistency 
with chiral perturbation theory, which describes the amplitudes at low energies well. 

In order to look into the characteristics of resonances, we need to investigate the partial- 
wave amplitudes. In order to compare the present results with the experimental data, we 
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normalize the transition amplitudes as follows: 

'--n^s),,. (18) 



Fig. 3 represents the partial-wave amplitudes for vrA^ scattering in the Sn and S31 channels. 
We find that the results with the MF and DF are similar to those with the DIM. The kink 
around 1500 MeV implies the resonance of A^*(1535). Hence, A^*(1535) is dynamically gener- 
ated in the form-factor regularization as well as in the dimensional one, once the subtraction 
parameters Oj are properly chosen. However, the present calculation fails to reproduce the 
higher resonance A^*(1650). Concerning A^*(1650), we mention that Ref. [l6i] has repro- 
duced it by introducing more parameters than the present work. Thus, in the present frame 
work, using the lowest-order Chiral Lagrangian with the on-mass-shell approximation is not 
enough to generate the A^*(1650) resonance dynamically. 

In Figl^ we draw the total cross sections of as functions of the laboratory momentum of 
71^ . Both regularization schemes describe them qualitatively well. However, there is some 
difference in detail, which is due to the fact that the corresponding propagators shown in 
Fig. n have different slopes as the energy increases. In particular, the total cross section of 
the TT^p iC^S" shows a noticeable dependence on regularizations. 

Now, we are in a position to look for the poles corresponding to the resonances in each 
channel of the S = 0, / = 1/2 meson-baryon sector. In order to find the pole corresponding 
to the resonances, we need to extend the T matrix to the complex energy plane in which 
the propgator has a cut starting from the threshold energy. 

We find that the pole is located at (1516 — 37i) MeV in the case of the DIM, while those 
with the DF and MF are at (1502 - Ui) MeV and at (1517 - Ali) MeV, respectively. These 



where two pions channel and p meson 



results are qualitatively similar to those in Ref. 
exchange are also included. 

We can express an approximated form of the T-matrix amplitude near the pole as follows: 

T,, ^ (19) 

where the residue Qi denotes the coupling strength to the resonance in the ith channel. 
Hence, we can easily determine the coupling strength of the resonance, using Eq. lfTIH) . The 
coupling strengths of A^*(1535) to each channel of S = 0,/ = 1/2 are listed in Table ITTTl We 



8 



find the following tendency of the coupling strengths: 



\9ttn\ < \9ka\ ~ \9vn\ < Is'keI 



(20) 



for all regularizations. With these results we can conclude that the A^*(1535) resonance is 



strongly coup 
those in Ref. 



ed to the KT, channel 



3. 



It is interesting to com par e the above results with 



14| . we get larger value of 



j. While almost all channels are similar to Ref. 
\gi\ for the KT, channel as quoted in Ref. [2^. The reason lies in the fact that the vrvrA^ 
channel may bring down the coupling strength of the KT, channel. 



IV. NUMERICAL RESULTS OF S = -1 MESON-BARYON SECTOR 

In this section, we investigate the S = — 1 sector, emphasizing on the resonances A(1405) 
(/ = 0) and A(1670) (J_= 0) in S wave. The A(1405) and A(1670) resonances have already 
been studied in Refs. jll |l5| in the chiral unitary model in which it was shown that A(1405) 
depends weakly on the choice of subtraction parameters while A (1670) is very sensitive to 
them . The fitting procedure is similar to that in Section IIIII We use /i = 630 MeV 
for the DIM and A = 1000 MeV for the DF, and MF. The average meson decay constant 
/ = 107.18 MeV for S = — 1 is adopted for all regularizations as in Ref. The fitted 
subtraction parameters are given in Table llVt where those of the DIM are taken from 
Ref. !^|. 

In the S = — 1, / = sector it is of great importance to look into the vrE mass distribution 
around the threshold of KN (~ 1432 MeV), since the A(1405) resonance arises from it. Fig. El 
shows the vrE mass distribution corresponding to the A(1405) resonance. As shown in Fig. El 
the peak is well reproduced as compared to the experimental data jl^. There is almost no 
difference between regularizations. 

The cross sections of the S = —1, = are plotted in FiglHl as functions of the lab 
momentum for K~p scattering. The results are in a fairly good agreement with experimental 
data. While difference between regularizations is found in the K~p K~p, K~p K^n, 
and K~p tt^S^ processes, there seems no dependence on the regularization in the case 
of the K^p — > K^n and K~p vr^A processes. 

We find three different poles in the scattering amplitude in the S* = — 1, / = channel. In 
particular the two poles around ~ 1400 MeV have been confirmed in recent studies based on 
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the chiral unitary model 
of the poles and coupling strength as done in Section IIIII Three different regularizations 
give almost the same results. 

In the S = —1, I = channel, the tendency of the coupling strengths are as follows: 

\9ke\ < \9va\ < \9kn\ < Is'tteI for pole 1, 
\9ke\ < \9riA\ ~ Ifi-TTsI < \9kn\ for pole 2, 
Is'tteI < Ifi-A'Tvl < \9r,A\ < \9ke\ for pole 3. 

V. DISCUSSTION AND SUMMARY 

We have investigated 5* = and S = —1 meson-baryon scattering, in particular, focusing 
on their dependence on regularizations. Starting from the effective chiral Lagrangian to 
the lowest order of the chiral expansion C\.jq, also known as the Weinberg- Tomozawa term, 
we have constructed the pseudo-potential which is used as a kernel of the Bethe-Salpeter 
equation. In order to solve the Bethe-Salpeter equation, we utilize the on-mass-shell ap- 
proximation. While the dimensional regularization is widely adopted in almost all works, 
we introduced three different schemes of the regularization: The dimensional regularization, 
the monopole type form factor (MF), and the dipole type form factor (DF). 

We examined first the dependence of the propagators on regularizations. While there is 
basically no difference among the regularizations in low-energy regions (below 1500 MeV), we 
found that in higher-energy regions (above 1500 MeV) the propagator with the dimensional 
regularization differs substantially from that with the monopole and dipole-type form factors. 
We fixed the subtraction parameters for the MF and DF at the threshold point. In addition, 
we used empirical values of the branching ratio and the partial wave amplitudes to fit the 
S = —1 and 5* = subtraction parameters for the form-factor regularizations. For the form 
factor regularizations, we have found that the on-mass-shell approximation is limited to the 
energy region lower that a certain energy due to the appearance of an unphysical branch cut. 
Therefore further study with the BS scattering equation solved numerically is necessary. 

We first calculated the partial-wave amplitudes for the S = meson-baryon sector for 
1 = 1/2 and I = 3/2. Though there exists difference between regularizations in higher- 
energy regions, the results were almost the same below 1500 MeV. The reason lies in the fact 



31, 



32, 



33 



3. 



In Table M and Table HVl we list the positions 
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that the propagator in the Bethe-Salpeter equation is rather sensitive to the regularization 
schemes in higher-energy region. We also calculated the total cross sections for the S = 
meson-baryon sector. While all three regularizatoins describe well them qualitatively, we 
found noticeable difference between regularizations in the total cross section of n^p i^'^S^. 
The position of the resonance A^*(1535) is similar to each other in all regularizaton schemes. 
However, the coupling strengths of A^*(1535) is rather sensitive to regularizations. 

The resonance A(1405) in the S = — 1 meson-baryon sector was well reproduced in the 
present work and turned out to be rather insensitive to the types of regularizations apart 
from small difference in higher-energy region, we also investigated the analytic structure of 
the partial-wave amplitudes and found that three poles exist in them. 

The regularization dependence of the the total cross sections become different according 
to the channels. While the processes K~p K~p, K~p K^n, and K~p 7rS+ depend 
on the regularization schemes, even in lower momentum region, those of K'p tt^'E^ and 
K~p — > p — > tt'^'E" show almost no dependence on regularizations. The process K~p — > 7r°A 
is changed by regularizatons only in higher momentum region. The coupling strengths 
strongly depend on regularizatons as in the S = — 1, / = meson-baryon sector. 

In the present work, we examined the dependence of meson baryon scattering on regu- 
larization schemes. We employed three different regularizations: The dimensional regular- 
ization, the form factor regularizations with monopole and dipole types. The differences 
due to regularizatons found in observables are mainly due to the different behavior of the 
propagator G{^/s) , because of which in higher energy region the regularizations change the 
prediction of the observables. Thus, we conclude that one need to vindicate methods used 
so far to describe meson-baryon processes. 
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FIG. 1: The real parts of the meson-baryon loop integrals for ttN {S = —1, I = 0) and ttS {S = 0, 

7 = 1/2) channels as functions of the center of mass energy Ecm(= \/s) for DF (solid), MF (dotted) 
and DIM (dashed). 
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FIG. 2: X 10-^2 as functions of E^m for KN (sohd), vrE (dotted), rjN (dashed) and KE (dot- 
dashed) channels. 



TABLE I: S = channels 





12 3 4 


7 = 1/2 


ttN rjN KA KT, 


7 = 3/2 


ttAT 7CS - - 
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TABLE II: Subtraction parameters for iS = meson-baryon sector for the DIM, DF and MF 





vriV r?iV KK 




DIM 


2.0 0.2 1.6 


-2.8 


DF 


3.00 0.93 2.38 


-1.87 


MF 


3.83 1.73 3.26 


-1.12 



1.0 




^ ■ ■ ■ ■ ■ ■ . — U.O ' ■ ■ ■ ■ ■ ■ ■ ' 

1000.0 1200.0 1400.0 1600.0 1800.0 1000.0 1200.0 1400.0 1600.0 laoo.o 



E cm [MeV] E cm [MeV] 

FIG. 3: 'kN scattering amplitudes {S = 0) for 6*11 and 5*31 channels as functions of the center of 
mass energy. Thin and Thick curves arc for real and imaginary parts, where the solid, dotted and 
dashed ones for DIM, DF and MF, respectively. Experimental data are shown by diamonds for 
the real parts and by circles for the imaginary parts. 



TABLE III: Coupling strengths l^^l of AT* (1535) to four different channels. 





■kN r]N KA KT. 


DIM 


0.46 1.79 1.09 5.17 


DF 


0.84 1.93 1.60 3.23 


MF 


1.00 2.02 1.78 3.54 
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FIG. 4: S = 0, /s = — 1/2 cross sections with different regularizations for ir^p K^T,'^, K'^A, and 
rjn as functions of the lab momentum Piab- The soUd curves are for the DF, the dotted one for the 
MF, and the dashed one for the DIM, respectively. 



TABLE IV: Subtraction parameters a for the S = —1 meson-baryon sector for the DIM, DF and 
MF. 





^KN 




«7rS 


0'r)K 






DIM 


-1.84 


-1.83 


-2.00 


-2.25 


-2.38 


-2.67 


DF 


0.22 


0.59 


0.47 


-0.12 


-0.20 


-0.36 


MF 


1.03 


1.36 


1.22 


0.64 


0.54 


0.36 
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FIG. 5: ttS mass distribution around the A(1405) resonance. The sohd curve represents the case 
of the DM, while the dashed and dot-dashed ones correspond to the case of the DIM and MF, 
respectively. Experimental data is taken from Ref. [2^. 

TABLE V: Poles of the scattering amplitude in the S = —1 and / = channel. 





Pole 1 (MeV) Pole 2 (MeV) Pole 3 (MeV) 


DIM 


1398 - 74z 


1429 - 14i 


1688 - 22i 


DF 


1392 - 7Ai 


1422 - 18i 


1663 - 9i 


MF 


1394 - 74i 


1425 - m 


1671 - 25i 



TABLE VI: Coupling strengths {gal for three poles. 







KN vrS ??A KE 




Pole 1 


1.43 2.06 0.54 0.45 


DIM 


Pole 2 


2.52 1.30 1.32 0.29 




Pole 3 


0.68 0.16 0.98 3.14 




Pole 1 


2.00 2.54 0.78 0.62 


DF 


Pole 2 


3.64 1.83 1.82 0.53 




Pole 3 


0.45 0.18 0.93 2.23 




Pole 1 


2.65 3.57 1.01 0.83 


MF 


Pole 2 


3.39 1.80 1.73 0.48 




Pole 3 


1.17 0.36 1.53 4.87 
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FIG. 6: S = —1, Is = cross sections with different regularizations for K~p —> K~p , K^n, tt^A, 

7r°E°, 7r+S~, 7r~S+ as functions of the lab momentum i^ab- The sohd curves are for the DF, the 

dotted one for the MF, and the dashed one for the DIM, respectively. 
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